Correspondence between HBT radii and the emission zone in non-central heavy ion collisions. 
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In non-central collisions between ultra-relativistic heavy ions, the freeze-out distribution is anisotropic, and 
its major longitudinal axis may be tilted away from the beam direction. The shape and orientation of this 
distribution are particularly interesting, as they provide a snapshot of the evolving source and reflect the space- 
time aspect of anisotropic flow. Experimentally, this information is extracted by measuring pion HBT radii as a 
function of angle with respect to the reaction plane. Existing formulae relating the oscillations of the radii and 
the freezeout anisotropy are in principle only valid for Gaussian sources with no collective flow. With a realistic 
transport model of the collision, which generates flow and non-Gaussian sources, we find that these formulae 
approximately reflect the anisotropy of the freezeout distribution. 
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PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq 



I. INTRODUCTION 

Femtoscopic intensity interferometry measurements use 
two-particle correlation functions to probe the space-time sub- 
structure of the system generated in heavy ion collisions, at 
the femtometer scale. At relativistic collision energies, exten- 
sive systematic studies have mapped out this substructure as a 
function of transverse momentum (pr), rapidity (y), collision 
energy an d me mass of the two correlated particles; 

see U for a recent review. These measurements probe the 
space-time geometry of the freezeout distribution- the distri- 
bution of last-scattering points of the particles. In addition 
to the source size and lifetime at freezeout, the momentum- 
dependence of the femtoscopic scales reveal the coordinate- 
space aspects of collective motion- "flow." 

It has long been recognized that measurements relative to 
the direction of the impact parameter of the collision are more 
sensitive to important underlying physics of the system, than 
are angle-integrated measurements. The azimuthal depen- 
dence of particle yields and spectra- often called "directed" 
and "elliptical flow"- are extensively used to extract the QCD 
equation of state (EoS) and transport coefficients of the quark- 
gluon plasma such as viscosity 12 0. On the other hand, 
azimuthally-integrated pj spectra can flag the existence of 
collective behavior, but are not as discriminating between dif- 
ferent models producing such behavior [e.g. HIE]. Similarly, 
the azimuthal dependence of jet quenching is a more dis- 
criminating probe of partonic energy loss, than azimuthally- 
integrated measurements |6]-[8]. 

It is possible that azimuthally-differential analysis might 
yield a similar improvement in sensitivity of femtoscopy. 
Only a few measurements of femtoscopic pion correlations 
relative to the impact parameter have been reported II91 4TT1 . 
though an extensive energy-dependence of these measure- 
ments is underway at RHIC fl2l . 

Most femtoscopic analyses of pion correlations extract so- 
called "HBT radii" (c.f. discussion in (I)) which fully de- 
scribe the emission distribution only if it is purely Gaussian 
and features no collective motion. Strictly speaking, neither 
of these conditions characterize heavy ion collisions, and there 



has been considerable activity in measuring the non-Gaussian 
features of the source via "imaging" techniques lfT3l [T4l that 
fit the source with a sum of spline functions. Large non- 
Gaussian tails are mostly explained by long-lived resonance 
production lfl"5l . 

Gaussian HBT radii seek to capture the bulk scales of 
particular interest; several studies exist, testing the corre- 
spondence between HBT radii and the source scales of non- 
Gaussian, flowing distributions from cascade calculations lfl6l 
im and blast-wave and hydrodynamical models [18]. 

Moving beyond HBT radii themselves, which characterize 
the geometry only of subsets ("homogeneity regions" lfl9l ) of 
the overall source, one is interested in the shape and orien- 
tation of the emission region as a whole. Retiere and Lisa 
have propsed l20l a formula connecting the azimuthal oscil- 
lations of HBT radii with the transverse anisotropic shape of 
the source. Here, we follow the same line and propose a for- 
mula for the tilt of the source relative to the beam direction. 
Both of these formulae are strictly valid only for Gaussian, 
non-flowing sources. In this paper, we test these formulae 
in the context of a realistic transport model featuring non- 
Gaussian freezeout distributions and strong flow. We find 
reasonable consistency between the source eccentricity and 
tilt as extracted directly from the space-time freezeout coor- 
dinates, and the same quantities estimated with the formulae. 
The discrepancy between the two provides an estimate of the 
systematic uncertainty one expects, when using azimuthally- 
differential pion correlation measurements to extract the un- 
derlying source shape and orientation. 

In section |nj we define the formalism and the connection 
between HBT radii and an anisotropic but simplified static 
Gaussian source. There, we present formulae connecting 
measurable quantities to the interesting features of the ge- 
ometry. In section III a realistic transport model (UrQMD) 



is used to generate a freezeout distribution featuring non- 
Gaussian geometry and strong collective flow. We discuss 
the non-trivial anisotropics of the distribution and the rela- 
tionship between regions of homogeneity |[T9l ED and the 
"whole" source of emission points. In section [TV] we sim- 
ulate an experimental analysis, using the UrQMD-generated 
distributions. We build two-pion correlation functions and 
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use the formulae presented in section [H] to estimate the 
source anisotropies. The calculations are compared to re- 
ported experimental results from Au+Au collisions at y/s/w = 
3.84 GeV, and predictions are given for heavy ion collisions 
at y/SNN = 30 GeV, relevant for FAIR and the RHIC energy 
scan. In section [V] we summarize. 



II. HBT RADII AND THEIR CONNECTION TO THE 
UNDERLYING SOURCE 

Femtoscopic two-particle correlation functions as a func- 
tion of the relative momentum q = p\ — p2 are often fitted in 
terms of Gaussian HBT radii R 2 , 




FIG. 1: (color online) The simplified parameterization of the freeze- 
out distribution in a heavy ion collision. In addition to the timescale 
and three spatial length scales, the ellipsoid may be tilted in the di- 
rection of the impact parameter x, relative to the beam axis z- 
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acterize the source, and the HBT radii are given by 



Indices i and j indicate the components of the q vector in 
the so-called Bertsch-Pratt "out-side-long" coordinate sys- 
tem [22-24 1 . Here, "out" points along the direction of the pair 
transverse momentum, and "long" points along the direction 
of motion of one of the incoming nuclei, say, in the direction 
of Au ions in the yellow ring of RHIC. We shall call this the 
"yellow nucleus" and its colliding partner the "blue nucleus." 
(We return to this arbitrary designation soon.) The "side" di- 
rection is given by the cross-product of "out" with "long." It 
is worthwhile to point out that, while one may simultaneously 
flip the signs of all components of q by swapping the desig- 
nation of particles 1 and 2, the correlation function depends 
only on even-order products of q's components; these prod- 
ucts have meaningful sign. 

We begin by considering a simple source of midrapidity pi- 
ons which is a Gaussian ellipsoid in space and time, with the 
major axis of the ellipse tilted with respect to the beam direc- 
tion, as sketched in Figure [T] The distribution is characterized 
by five parameters: a temporal scale, three spatial scales and 
a tilt angle 
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where the primes on Gy and a z i signify that these correspond 
to principle axes of the ellipse. 

Its transverse eccentricity about its (tilted) major axis is de- 
fined as 
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where px and P/ are the transverse and longitudinal velocities 
of the pion pair. "Cross-term" radii R 2 , ■ vanish by symme- 
try (25] |26). 

For the more general case, there are six HBT radii, and they 
depend on the azimuthal angle (j). This angle is meaningfully 
defined over the range [0,27t] about the beam direction relative 
the impact parameter, which is defined as the direction per- 
pendicular to the beam, pointing from the yellow nucleus to 
the blue one. Swapping the designations of the "yellow" and 
"blue" nuclei reverses the direction of the impact parameter (x 
direction); however, it also reverses the "long" (z) direction, 
defined earlier. Hence, the sense of the tilt 0s is well-defined; 
a source with positive tilt features a positive spatial x — z cor- 
relation, as shown in figure [T] Experimental measurement of 
the sense of the tilt [9| requires measuring the direction of the 
impact parameter, for example through directed flow of net 
baryons at forward rapidity. 

The HBT radii, measured as a function of angle (j) relative 
to the beam axis, are driven by source widths o v , o Y , a z , rather 
than a x i,G y ,O z i, The relationships between these widths are 
given by 



a 2 = a 2 cos 2 0, s + a 2 sin 2 0, s + a 2 , sin 2Q S 
a 2 = a 2 sin 2 S + a 2 cos 2 Q s - a 2 , sin 20 s 



(5) 



If S = and a x = o y = a± (or equivalently in an 
azimuthally-integrated analysis), only three parameters char- 



where a 2 , is the covariance between x and z in the source func- 
tion; see [25] for details. 
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The HBT radii are related to these widths as 
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In analogy with Equation [3] we identify the eccentricity of 
the source around the beam axis as 



If the first-order plane is identified, first-order azimuthal os- 
cillations in R 2 sj and R 2 ol are measureable. In this case one ob 
tains the tilt angle according to 1 25 1 

0, ; tan 1 ( ... I;- (II) 
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The transverse eccentricity in the "natural" frame tilted rela- 
tive to the beam axis is 

2R 2 s2 (1 + cos 2 9. s ) + (r 2 - R 2 Q ) sin 2 Q s - 2R 2 n sin 29* 



R 2 (1 + cos 2 9. s ) + [2R 2 2 +R 2 j sin 2 9, + 2R 2 sll sin 29, 

(12) 
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Equations [6] express the explicit 0-dependence of the HBT 
radii for a non-flowing source; all other variables in the equa- 
tions are constants, for a non-flowing source. For a source 
with flow, the constants, e.g. o y , may themselves depend on 
0. In this case, our equations will be violated; below, we quan- 
tify this violation and its effect on the extraction of 9 S and e. 

Experimentally, one measures the squared HBT radii and 
calculates the Fourier coefficients quantifying their azimuthal 
dependence, per 
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In our no-flow Gaussian model, then, the source geome- 
try and orientation are extracted from the Fourier coefficents. 
Identifying the tilt angle requires [25 26 1 measuring HBT rel- 
ative to the first-order reaction plane Q. Published results 
from the STAR OS) and CERES Ol collaborations use only 
the 2 nd -order plane and so report the source eccentricity only 
about the beam axis. (Efforts to perform the analysis relative 
to the first-order plane are underway at RHIC |[T2l .) In this 
case S2Q1 
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For the moment, we ignore the tildes above the n ^ 1 Fourier 
coefficients here and below. They represent a trivial finite- 



Ill. FREEZEOUT DISTRIBUTIONS FROM A MORE 
REALISTIC MODEL 

Our simple model of the emission function given in equa- 
tion|2]is unrealistic in at least two ways. Firstly, while realistic 
emission functions are often roughly Gaussian, they are never 
perfectly so; in this case, the two-pion correlation function 
is likewise non-Gaussian. Extracting Gaussian HBT errors 
through fits with equation [T[ then, could in principle generate 
considerable mis-representation of the emission function. 

The second over-simplification of the source discussed 
above is its lack of collective flow, which generates correla- 
tions between a particle's emission position and its momen- 
tum lfl9l . For example, an explosively flowing source will 
boost particles emitted from its right side, towards the right. 
The freezeout distribution of particles with a given momen- 
tum vector is known as the region of homogeneity for that 
momentum vector. In principle, the homogeity regions for dif- 
ferent azimuthal angles might be completely disjoint and un- 
related, obviously invalidating equations |6} [T0| [TT| and [T2"| In 
practice, the homogeneity regions in blast-wave models [20] 
or hydodynamic simulations [27 1 are naturally related. For 



binning correction discussed later in section IV 



boost-invariant sources, equation 10 is estimated to be good 
to ~ 30% in these models [20|. 

While blast-wave and boost-invariant hydro models do fea- 
ture non-Gaussian sources and collective flow, they are still 
simplistic. Firstly, any boost-invariant model by definition has 
no tilt relative to the beam axis; thus they are unable to ac- 
cess physics associated with 9,. Secondly, they typically use 
the Cooper-Frye formula to model particle freezeout from a 
calculated or parameterized hypersurface; while momentum- 
space observables (e.g. V2) may be insensitive to this pro- 
cedure, interferometry is known to be quite sensitive to the 
freezeout proceedure. 

In this section, we use the Ultra-Relativistic Quantum 
Molecular Dynamics Model (UrQMD 3.3)||46) to generate a 
realistic freeze-out distribution with fully three-dimensional 
dynamics [28, 29 1. UrQMD is a covariant transport approach 
to simulate the interactions between hadrons and nuclei up to 
relativistic energies. It is based on the propagation of nucleons 
and mesons accompanied by string degrees of freedom with 
interaction probabilities according to measured and calculated 
cross sections for the elementary reactions. Hard scatterings 
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FIG. 2: (color online) Pion emission points from UrQMD simula- 
tions of Au+Au collisions with collision energy v/*W = 3.84 GeV 
(corresponding to a 6 AGeV beam incident on a fixed target) and 
impact parameter b = 4 — 8 fm. Colored contours (identical for the 
three panels) show the emission point density of all pions with pj < 
0.4 GeV/c. Black contours in the upper, middle and lower panels 
indicate the density of emission points for pions with (]) = [— I , 5 1 , 



[f ,f ] and [f,^], respectively. 



with large momentum transfer are treated via PYTHIA. For 
detailed comparisons of this version to experimental data, 
the reader is referred to |30|. Previous HBT studies with 
UrQMD have been reported in Il3lti34l . 

In UrQMD, the freeze-out space-time point is naturally de- 
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FIG. 3: Freezeout distribution of pions from y/Jfm = 3.84 GeV 
(6 AGeV beam energy on fixed target) Au+Au collisions with im- 
pact parameter b = 4 — 8 fm in the reaction (x — z) plane, as calculated 
from UrQMD. 



fined as the last hard interaction of a particle. The freezeout 
distribution of pions from y/sim = 3.84 GeV Au+Au colli- 
sions with impact parameter b = 4 — 8 fm in the reaction plane 
{x — z plane) is shown in figure [3] The source has an obvious 
tilt structure relative to the beam axis. We use the parameteri- 
zation of equation[2]to fit the three-dimensional freezeout dis- 
tribution for all pions with pj < 400 MeV/c and \y\ < 0.6- not 
only those at a given angle (j). From this direct analysis of the 
freezeout coordinates- obviously not possible in experiment- 
we find the parameters listed in the third column of table [I] It 
is clear from figure[3]that x — z correlation in the freezeout dis- 
tribution has structure that cannot be captured in a single tilt 
number; indeed, the tilt is scale-dependent, growing as one 
focuses on the peak of the distribution. This "twist" feature, 
which has been observed in simulations before [25 1, might be 
physically interesting and experimentally accessible; we leave 
exploration of this effect for a future work. For the purpose 
of this work, we identify a range of tilts arising from fitting 
equation [2] to the distribution and varying the fit range in co- 
ordinate space from 10 fm < Ax, Ay, Az < 40 fm. This leads to 
the range shown in the left column of table [I] 

Dynamics naturally leads to a strong correlation between 
a particle's final momentum and the freezeout position; ho- 
mogeneity regions are naturally reflected in the final state. 
Figure [2] shows homogeneity regions from UrQMD in the 
x — y plane for particles emitted at three azimtuthal angles. 
The homogeneity region for particles emitted at a given angle 
(|)i clearly differs from that for particles emitted at §2 ^ <j>i. 
Thus, in addition to the explicit (^-dependence of the HBT 
radii seen in equations [6] there is an additional implicit de- 
pendence [35—38 1. HBT radii measured at a given momentum 
vector (§,pr,y) probe only the geometry of the homogeneity 
region for that momentum vector. A priori, it is far from clear 



that equations 10 1 1 and 12 which attempt to relate HBT ra- 
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dius oscillations to the geometry of the "whole" source, will 
prove good approximations. 

In the next section, we test these relations- true for the 
simplest toy model of Equation [2]- with UrQMD-generated 
freezeout distributions. 



IV. HBT RADII AND ANISOTROPY PARAMETERS FROM 
URQMD-GENERATED CORRELATION FUNCTIONS 

We start by generating two-pion correlation func- 
tions, analogous to those measured experimentally, from 
UrQMD events. We will then proceed to fit these correla- 
tion functions with the Gaussian anzatz of equation [T] as is 
done in experimental analysis. Finally, we extract Fourier co- 
efficients characterizing the oscillations of the HBT radii with 
angle; from these we extract the source anisotropy parameters 
that would be obtained in an experiment. 

To simulate experimental conditions, two-pion correlation 
functions were constructed from the UrQMD events using 
the so-called weighting method l47l . In this method, pairs 
of identical pions are selected according to a Monte Carlo al- 
gorithm; the correlation function in a given (q r) ,q s ,qi) bin is 
equal to the pair-wise average of the squared two-pion squared 
wavefunction. Considering only quantum symmetrization ef- 
fects, the correlation function is computed as 

C(q ,q s ,qr,$) = (l+cos(-^Ar iM )> ( | ) , (13) 

where q = p\ — pi is the relative momentum and Ar is the 
space-time separation of the particles at freezeout. 

As explicitly denoted in equation [13] the correlation func- 
tions were generated for 8 45°-wide bins in pair angle (]) = 
Z(KT,t>), where Kj = (pr.i +Pt,i) /2 is the average trans- 
verse momentum vector of the pair. Hence, for — | < <|) < | 
(4r < (]) < ^), pions from sub-region indicated by the black 
contours in the top (bottom) panel of figure[2]are used to con- 
struct the correlation function. 

In figure |4] are plotted, for each bin in (j), two-dimensional 
slices of the three-dimensional correlation function in the 
q„ — q s , q s — qi and qi — q H planes; in each case the unplotted 
relative momentum component q\ < 4 MeV/c. For a represen- 
tative (j) bin, one-dimensional slices of the correlation func- 
tion in the out, side, and long directions are shown in Fig- 
ure [5] Most femtoscopic correlation analyses focus on the 
one-dimensional projections, since the correlation in three- 
dimensional space factorizes; that is, there is no covariance 
between components qi and qj^ in the correlation function. 
This is not the case when the analysis is performed differen- 
tially in (j) |9, 25, 26, 35-38), as is clear from the tilted struc- 
tures in g-space seen in figure |4] These tilts in the individual 
correlation functions in g-space for a given angular selction 
in K are not to be confused with the overall spatial tilt of the 
source sketched in figure [T] 

As in an experimental analysis, the correlation functions are 
fitted with the Gaussian functional form of equation [T] Two- 
and one-dimensional slices of these fits are superimposed on 
the correlation functions in figures [4] and [5] The six resulting 
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FIG. 4: (color online) Pion correlation functions from Au+Au col- 
lisions with collision energy v/*MV = 3.84 GeV and (corresponding 
to a 6 AGeV beam incident on a fixed target) and impact parameter 
6 = 4 — 8 fm, as calculated with the UrQMD model. Projections 
in the q — q s (left column), q a — q\ (middle) and q s — qi spaces are 
made, with the unplotted ^-component smaller than 4 MeV/c. Corre- 
lation slices are shown for pion pairs in 45°-wide bins centered at 
angles indicated to the right of each row. Shaded (color online) con- 
tours represent the calculated correlation function. Black contours 
represent 2-dimensional slices of the three-dimensional Gaussian fit 
to the correlation functions for each selection in <|). 



HBT radii are plotted as a function of (j) on figure [6] As in a 
standard azimuthally-integraged analysis, the "diagonal" radii 
Rj, j = o,s,l are driven by the width of the correlation func- 
tion in the direction j. The "cross-term" radii Rfiuf quantify 
the correlations between q components- the tilts of the corre- 
lation function; e.g. the ^-dependence of the tilt of the corre- 
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FIG. 5: One dimensional slices in the three components of relative 
momentum, for pions emitted at \§\ < 22.5°. Curves represent one- 
dimensional slices of the three-dimensional Gaussian fit to the corre- 
lation function. 




FIG. 6: The six HBT radii extracted from Gaussian fits (equation[TJ 
to the correlation functions for eight selections on (]). Curves rep- 
resent a Fourier decomposition (equation [5), including terms up to 
second order, where the Fourier components are determed according 
to equation|9] See text for details. 



lation function in the q s — q\ space seen in the right column of 
figure pleads directly to the first-order oscillation of R 2 sl seen 
in figure [6] Figure [7] shows the radii for ^Js^n =7.7 GeV 
collisions. 

The curves on figures [6] and [7] represent equations [9] using 
Fourier coefficients R^ n calculated according to|9] 

Before we apply equations T0p2 we note that the argu- 
mentation of section|n implicitly assumed that the R 2 (())) were 
measured as a continuous function. In reality, in our analysis 
as in experiment, the correlation function is measured for (]) 
within bins of width A(j). The amplitude of the « th -order os- 
cillation of a binned function is reduced from that of the un- 
derlying function. To correct for this finite-binning artifact, 
we calculate the underlying ("true") Fourier coefficients f& n 
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FIG. 7: Same as figure|6] but for collisions at ^/snn = 7.7 GeV. 



from the ones extracted from the binned radii by 



Rl, 



nA<|>/2 
sin 



Rt. 



(14) 



= 1.026^ ^ and ^ 2 = 1.11 12^ 2 . 



For our 45° -wide bins 

These binning-corrected Fourier coefficients are used in equa 
tions [T0l[T2l ED 

As in most experimental analyses, the correlation functions 
from UrQMD simulations are not purely Gaussian, since the 
source itself is not Gaussian in coordinate space, due to space- 
momentum correlations (flow), resonance contributions, etc. 
Of special interest for the present study is the additional fact 
that the source is not characterized by a tilt angle independent 
of spatial scale- the "twist" discussed in section [III] Follow- 
ing a standard experimental approach [e.g. 39], we perform a 
fit-range study, in which we vary the range in q () , q s , qi, over 
which we perform the fit with equation [T] In particular, we fit 
the correlation functions in the range —q ma x < qo,qs,qi < qmax 
for q max = 60 — 150 MeV/c. The resulting anisotropy param- 
eters for y^SNN = 3.84 GeV collisions, calculated according 
to equations [T0p2] are shown on figure [8] The dependence 
of 05 on q mux is readily understood. Large values of q probe 
smaller values of coordinate space; thus, as q max is increased, 
the fit is increasingly sensitive to the large tilt structure seen at 
small scales in figure [3] Figure [8] then, is itself a measure of 
the "twist" structure in coordinate space, though there may be 
more sophisticated ones. For our purposes, however, we take 
the variation of the anisotropy parameters plotted in figure [8] 
to define a range of values one might expect from an experi- 
mental study. Since a typical experimental analysis would fit 
at least out to 100 MeV/c (in order to include all of the peak 
signal), the value ranges listed in the fourth column of table [I] 
correspond to 100 MeV/c < q max < 150 MeV/c. 

The agreement with the parameters extracted via direct 
analysis of the UrQMD freezeout coordinates is fair; and we 
discuss this further in the next section. 

The fifth column in table [I] lists the anisotropy parame- 
ters based on HBT radii measured by the E895 collabora- 
tion in a fixed-target experiment at Ej, eam = 6 AGeV J5]. 
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FIG. 8: Source anisotropy parameters extracted from two-pion corre- 
lation functions for UrQMD-generated Au+Au collisions at = 
3.84 GeV at b = 4 — 8 fm, are plotted as a function of the range in 
relative momentum over which the correlation functions are fitted. 
Respectively, top, middle and bottom panels show the source tilt an- 
gle 05, eccentricity about the beam axis e and eccentricity about the 
tilted axis e'. See text for details. 
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UrQMD: HBT 
oscillations 
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3.84 GeV 
(6 AGeV) 
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0.13-0.17 
34° -41° 
0.21-0.26 


0.095-0.096 
27.4° -27.9° 
0.200 - 0.206 


0.30±0.15 

26° ± 7° 
0.38±0.19 


7.7 GeV 
(30 AGeV) 


£ 
£' 


0.11-0.14 

14° -21° 
0.12-0.16 


0.090-0.091 
U.7°-11.9° 
0.109 





The UrQMD calculation reproduces the tilt angle very well, 
the eccentricities somewhat less well, though experimental 
uncertainties are large. UrQMD calculations in the bot- 
tom row of the table represent a prediction for collisions at 
\/snn = 1.1 GeV, which will be measured at FAIR and RHIC. 



TABLE I: Source anisotropy parameters, for Au+Au collisions at two 
collision energies, with b = 4 — 8 fm for pions with pj < 0.6 GeV/c, 
and \y\ < 0.6. Impact parameter and momentum cuts were selected 
to match published data from the E895 collaboration |9|. Third col- 
umn: estimates from a Gaussian fit (equation [2) to the freezeout 
distribution from UrQMD events. Fourth column: estimates us- 
ing equations 1 1 0| 1 2| on the azimuthal oscilations of HBT radii from 
UrQMD events. Fifth column: same as column four, but using ex- 
perimental data from E895. Experimental data at 7.7 GeV will be 
analyzed at FAIR and RHIC. 



V. DISCUSSION AND SUMMARY 

The connection between anisotropic particle emission dis- 
tributions and femtoscopic correlation functions has been dis- 
cussed in detail. Gaussian fits to three-dimensional two-pion 
correlation functions result in six HBT radii (in contrast to the 
usual three), all of which depend on the pair emission angle 
relative to the impact parameter. The order and magnitude of 
the radius oscillations are quantified by Fourier coefficients, 
which may be directly related to the spatial tilt and eccentric- 
ity of a Gaussian source with no space-momentum correla- 
tions (arising, for example, from collective flow). For a more 
realistic freezeout distribution, the connection between the os- 
cillating radii and the source anisotropy is only approximate. 
Using a sophisticated transport model to simulate the entire 
collision evolution and particle decoupling, we have studied 
the freezeout distribution in coordinate space. The "whole" 
distribution indeed features a tilted structure and eccentricity, 
but shows other, less trivial anisotropic structures as well. 

The experimentalist, of course, cannot study the distribu- 
tion directly in coordinate space. We have simulated the 
experimental situation by constructing two-pion correlation 
functions for different pair emission angles, fitted each with 
a Gaussian functional form and extracted Fourier coefficients 
for each radius. These coefficients were then used in the for- 
malism derived in section[n]to estimate the tilt and eccentric- 
ity. The agreement between these estimates of the anisotropy 
and the more direct study in coordinate space, was fair. Pre- 
vious studies ll20l of boost-invariant hydrodynamic and blast- 
wave models found that the eccentricity values estimated from 
two-particle radii were within 30% of the "true" values; this 
has been used as a systematic error for the eccentricity in ex- 
perimental studies [ 10 1 . Our results using UrQMD are consis- 
tent with this 30% value. The present study is the first estimate 
of the corresponding uncertainty in the tilt, 9^; the agreement 
is somewhat worse, on the order of 35%. Given the compli- 
cations of dynamically-induced homogeneity regions, a time- 
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evolving emission distribution, non-Gaussianness and "twist" 
effects, one might easily have expected much worse agree- 
ment. 

However approximate, quantifying the connection between 
the radius oscillations and the underlying geometry can be 
useful. Ideally, the correct model of a heavy ion collision will 
reproduce all experimental observations; here, this means the 
HBT radii and their dependence on azimuthal angle. How- 
ever, when observations are reproduced and others not, con- 
nections to the underlying scenario are important. For exam- 
ple, if a model reproduces Ri ong and R s u e but over-estimates 
Rout [e-g- 40 -42 1, attention immediately turns to emission du- 
ration, which may be associated with the nature of the tran- 
sition between confined and deconfined states, latent heat, 
etc I22ll23ll43l. 

For the azimuthal dependence of HBT radii, the tilt an- 
gle and eccentricities probe different aspects of the dynam- 
ics. At AGS energies (~ 3.5 GeV), B$ reflects the dynamics 
behind directed flow |25| in the earliest stage of the collision 



and shows strong dependence of the equation of state used in 
transport calculations [9|. Meanwhile, the eccentricities rep- 
resent the geometric and temporal fl4ll aspects of elliptic flow. 
It is hoped that the connections we have discussed here, be- 
tween experimental observations and the underlying source 
anisotropy, will help future comparative studies focus on the 
physics most relevant to each observable. 
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